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A global analysis of the data from all the solar neutrino experiments combined with the recent 
KamLAND data is presented. A formula frequently used in the literature gives survival probability 
for three active solar neutrino flavors in terms of a suitably-modified two-flavor survival probability. 
Corrections to this formula, which depend on #13 and Sm^, are calculated. For the mass scale 
suggested by the atmospheric neutrino experiments the contributions of Sm^i to these corrections is 
found to be negligible. The role of #13 in solar neutrino physics is elaborated. For electron neutrino 
oscillations into another active flavor, we find best fit values of tan 2 612 ~ 0.46, tan 2 #13 ~ 0, and 
5m|i ~ 7.1 x 10~ 5 eV 2 . It is found that the combined solar neutrino and KamLAND data provide 
the limit cos 4 #13 > 0.8 at the 90 % confidence level. 
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I. INTRODUCTION 



Seminal developments in neutrino physics took place during the last few years. Observation of the charged-current 
solar neutrino flux at the Sudbury Neutrino Observatory (SNO) pj together with the measurements of the ^©-electron 
elastic scattering at the SuperKamiokande (SK) detector 2] established that there are at least two active flavors of 
neutrinos of solar origin reaching Earth. Furthermore this SNO charged-current, and subsequent neutral-current Q 
and day-night difference 3 measurements confirmed the prediction of the Standard Solar Model (SSM) for the total 
neutrino flux jaL Already in the summer of 2002 global analyses combining SNO and SuperKamiokande data with 
earlier chlorine @ and gallium 0, HQ measurements indicated the so-called large mixing angle (LMA) region of the 
neutrino parameter space as the most likely solution [To| . 

In December 2002 the KamLAND reactor neutrino experiment announced its first resu lts Ji ll . The KamLAND data 
also points to the LMA region, a result which was confirmed by independent analyses [H ElO E E H E H3 . 
One aim of our paper is to present a global analysis of all solar neutrino data along with the recent KamLAND data to 
examine neutrino parameter space. Our other goal is elaborate on the role of the mixing angle between the first and 
the third neutrino families, #13, in solar neutrino physics. Currently the value of #13 is probably the most-pressing open 
question in neutrino phy sics. Since the CP- violating phase appears together with #13, understanding CP violation in 
neutrino oscillations |2lJ, 1 2 J . |23| requires precise knowledge of 6*13 . A priori one expects the impact of #13 on solar 
neutrino physics to be minimal. However after recent significant accumulation of neutrino data we can now provide 
a more precise assessment of this conjecture. 

A formula frequently used in the literature gives survival probability for three active solar neutrino flavors in 
terms of a suitably-modified two-flavor survival probability. We calculate systematic corrections to this formula in 
Section II. These corrections depend on #13 and Sm^, however or the mass scale indicated by the atmospheric neutrino 
experiments the contributions of Sm^ to these corrections is found to be negligible. We use this formula to investigate 
the role of #13 in solar neutrino physics. Details of our statistical analysis are presented in Section III. Finally Section 
IV contains a discussion of our results. 



II. MATTER-ENHANCED NEUTRINO OSCILLATIONS 

We denote the neutrino mixing matrix by U a i where a denotes the flavor index and i denotes the mass index: 

* Q = XXi*i- (!) 
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For three neutrinos we take 



U a i — T23T13T12 = 



'10 

C23 S23 

,0 —523 C23 




(2) 



where C13, etc. is the short-hand notation for cos0i3, etc. Note that individual matrices, not their matrix elements, 
are called T23, T13, and T12, respectively in Eq. ■ I n ^13 the notation S^ 3 was used to indicate (sin0i3)e 1 ^ where 
4> is the CP-violating phase. We will ignore this phase in our discussion. The evolution of the three neutrino species 
in matter is governed by the MSW equation [24], I2M I2T 
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where the Wolfenstein potentials in neutral, unpolarized matter are given by 

V c (x) = V2G F N e (x) 

for the charged-current and 




(3) 



(4) 



V n {x) 



V2 



G F N n {x). 



(5) 



for the neutral current. Since V n only contributes an overall phase to the neutrino evolution in the Sun we will ignore 
it in the rest of t his p aper. 

Following Rcf. 27] we introduce the combinations 



IV = cos 6> 23 # M 
* r = sin 6» 23 IV 



- sin 2 3*r, 
COS 023 ^r, 



(0) 
(7) 



which corresponds to multiplying the neutrino column vector in Eq. with T^ 3 from the left. Eq. © now takes 
the form 




/E 1 ' 
T13T12 E 2 
V E 3 . 



7~>t 7~M 
1 12^13 





(8) 



where we used the the fact that T23 commutes with the second matrix on the right containing the matter potential. 
Note that the initial conditions on ^ M and ^> T are the same as those on VP^ and 'J,-: They are initially all zero. 
Consequently the mixing angle 023 does not enter in the source- and detector-averaged solutions of the matter- 
enhanced neutrino oscillations. 



We next perform the transformation 





COS0i 3 *e +sin0i3* T 
-sin 013^ + cos0i3* T 



after which Eq. (JSJ) takes the form 



(9) 
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'at 

where we dropped a term proportional to the identity. In this equation 7i is given by 

-±Ksin20 13 



\V- A 2 icos20 12 
H= I lA 2 i sin 20 12 
-kVrSin 20i 



|A 2 isin20 12 

+ A 2 1 COS 2012 
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where we introduced the modified matter potential 







i(A 31 + A 32 ) + K 



2 Y • 



(10) 



(11) 



V = V C cos 2 13 



(12) 
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and the quantity 



2E 



5m% 
2E 



(13) 



The corresponding mass matrix was considered in Ref . [23] • 

Eq. (|ll)fl describes the exact evolution of the neutrino eigenstates through matter. Assuming the mass hierarchy 
7713 > m,2 > mi, there are two MSW resonances. The lower-density resonance is apparent in this equation. Although 
it is not immediately obvious this equation also correctly describes the higher density resonance. One observes from 
Eq. © that the appropriate initial conditions are 



'<p e (t = 0)\ 

¥>/.(* = 0) = 

^r(* = 0)/ 



COS 6*13 


- sin 6*i3 , 



(14) 



These initial conditions need to be satisfied where the neutrino is produced, which is not necessarily at the center of 
the Sun. 

We look for the solution of Eq. (|1L)|) appropriate for the solar neutrino problem. Specifically we assume that sin #13 
is small as indicated by the reactor neutrino experiments prior to KamLAND |29l 1.3(1 131| . We also interpret the value 
8m\ 2 = 3x 1CU 3 eV 2 observed by the atmospheric neutrino observations at the SuperKamiokande detector j^] to 
imply that 5m\ 2 ~ ^ m ai since the LMA value of &m\ x deduced from the solar neutrino experiments is lower by more 
than one order of magnitude. Hence we take A31 ~ A32. Using the approximate electron density given by Bahcall 
331 



N e (r) = 245 exp(- 10.54r/R & )N A cm" 



(15) 



where Na is the Avogadro's number and i? Q is the radius of the Sun, we write the potential of Eq. (0J in the 
appropriate units as 



Using the value Sm^ 
condition 



V c {r) = 1.87 x 1(T 5 x exp(-lO.54r/i? ) eV 2 /MeV. (16) 
10~ 3 eV 2 and the Eq. Ijltjll . even for the highest energy (~ 15 MeV) solar neutrinos the 



Voir) 
a 3 i 



< 1 



(17) 



is satisfied everywhere in the Sun. Indeed as one moves away from the center of the Sun this ratio is much smaller 
especially for the lower-energy solar neutrinos. We will use this ratio as a perturbation parameter in solving the 
neutrino propagation equations. 

The differential equation involving the derivative of <p T is 



aip e + bip T 



where we introduced the abbreviated notation 



-V c sin 6»i3 cos 6»i3 



and 



1 



(A 3 i + A 32 ) + U c 



A31 - \V C {1 - 3 sin 2 13 ) = A31 - e. 



2 V 2 2 

Using the initial conditions of Eq. (|14|l . Eq. H18|) can be immediately solved to express <p T in terms of cp e 

r* 

T + i cos 0i3 



/ n . a —if b{t')dt' 

<p T {t) = sin 0i3 e J o 



dt'v c (t')cp e (tyIo b(t " )dt " 



(18) 



(19) 



(20) 



(21) 



Of course in general one still needs the value of ip e everywhere in the Sun. However when the limiting condition of 
Eq. I|17f) is satisfied one can use a technique employed in Ref. |34| to obtain an approximate expression. Writing 
de lA3lt = iA3ie lA3lt dt one can integrate the integral in Eq. 121fl by parts once to obtain 



■ f i-Lfjr i V b(t")dt" Vc{t — 0) 1 f jA„,t' ^ \tz r±r\ i V e(t")dt" 

i dt'V c (p e e Jo y ' = 1- ! C os0i3--— / e^ 3lt — V c (t')tp e (t')e Jo <■ > 

Jo A31 A31 J dt 1 



(22) 
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where we used the initial condition of Eq. I|14fl on ip e and the fact that V c is zero once the neutrinos leave the Sun. 
In the last integral in Eq. (IL'l'li we ignore the derivative of V c since it is many orders of magnitude smaller than V c 
itself at the core of the Sun (cf. Eq. ifTfljl ). We calculate the derivative of <p e using the Eq. l|Tn]l. Integrating this last 
integral by parts again we obtain an expansion of (p T outside the Sun: 



tp T (t) = sin 6»i3 e 

-1 + COS 2 013 



i £ Ht')dt' 



~t 1 



'31 



e (l-2sin 2 0i3 



where we defined 



V c (t = 0) 



V31 



(23) 



(24) 



Note that Eq. [12: il) does not represent <p T inside the Sun since we have already taken the limit V c = 0. 

The differential equation one can write from Eq. (|10|l involving the derivative of <p e includes a term containing tp T 
( — 7}V C sin 20i3 <Pt)- However substitution of ip T into this term using Eq. i|23|) indicates that this term is second order 
in sin 0i3 and can be discarded in a calculation which is first order in sin0i3. In this approximation the evolution of 
the first two components in Eq. I|1U[) decouple from the third yielding 



d 



dt \<Pn 



2 ' 



A 2 i cos 20i2 
|A 2 i sin 20i2 



|A 2 i sin 20i2 
-i^ + A 2 i cos 20i2 



(25) 



which is the standard MSW evolution equation for two neutrinos except that the electron density N e is replaced by 
iV e cos 2 013. The fact that one needs to solve the evolution equations with this new density was already emphasized 
in the literature [35ll3r| . Since this 2x2 "Hamiltonian" in Eq. (|1(J[) is an element of the SU(2) algebra, the resulting 
evolution operator is an element of the SU(2) group and the evolving states can be written as 



<p e (t)\ = ( $ e (t) 



Kit) 



COS 013 




= U 



COS 013 





where & e {t) and $ M (t) are solutions of Eq. (fTT)f) with the initial conditions $ e (i = 0) = 1 and ^(t = 0) = 0. 
We can write the electron neutrino amplitude using Eq. 10 as 

<f e = cos 0i3 ip e - sin 0i3 <p T . 

Substituting <p e from Eq. $Zffy and ip T from Eq. into Eq. (|77|) we obtain 



cos 0i3 $ e - sin 0i 3 e 



dt' 



D, 



(26) 



(27) 



(28) 



i f* b(t')dt' 
c : JO ; 



where we defined D = (p T e Jo " l ' /sin0i3 for convenience of notation, 
electron neutrino survival probability in our approximation: 

P 3x3(^e -» V e ) = COS 4 013 ^2x2(^6 ~> V e with N e COS 2 i3 ) 



Squaring Eq. (|28|l we finally obtain the 



sin 4 0i3 



l + 2ecos^0i 3 1 



cos 20i 



+ f COS 4 013 



, Sm 21 
' Sm 31 



COS 2012 — 1 



(29) 



f in this equation is given by Eq. (|24H . P2x2(^e 



with N e cos 2 013) is the standard 2-fiavor survival probability 



calculated with the modified electron density N e cos 613 and the standard initial conditions ($ e = 1 and <i> M = 0). 
Either a number of exact [39l I40L ItH \I$L or approximate 0, EBl |46| solutions of the neutrino evolution equations 
available in the literature or numerical methods can be used to calculate the survival probability of Eq. (|29|l . 

For £ = Eq. I|29l) reduces to the formula widely used in the literature (see e.g. Refs. 36, 38] ). We show the 
relative contributions of various terms to the survival probability given in Eq. 129|) in Fig. ^ We picked values of 
the neutrino parameters representative of the solar neutrino results in this figure. The upper solid line is the total 
survival probability numerically calculated for the 3-flavor mixing by solving the neutrino evolution equation, Eq. 
©, exactly using the method of Ref. 01 • The filled squares on top of the figure represent the contribution of 
the term proportional to cos 4 0i3 alone (the first term in the right side of Eq. (|29fl ) to the survival probability in 



5 



the approximation described above. One observes that this term alone is an excellent approximation to the total 
probability. To illustrate this we plot separately the contribution of the term proportional to sin 4 # 13 alone to the 
survival probability (the long-dashed line in the middle). Clearly the second term is approximately three orders of 
magnitude smaller than the first one. The correction terms proportional to the parameter £ are even smaller by two 
more orders of magnitude indicating the validity of our expansion in terms of the parameter £. It is worth emphasizing 
that the quantity Sm^ enters the survival probability in Eq. 1(29(1 only through the terms proportional to £, hence its 
effect is minimal: For solar neutrinos the 3-flavor survival probability depends only on two mixing angles, #12, #13, 
and one mass-difference squared, Sm^. Having gained confidence in the validity of Eq. 1)29(1 . we use it in our analysis. 



III. STATISTICAL ANALYSIS 



There is an extensive literature descri bing methods to calculate the goodness of a fit and confidence levels of 
allowed regions (See e.g. Refs. |4^ . l49ll5?l l5l| and other references we cite below). In our global analysis, we use the 
"covariance approach" , in which least-squares function for solar data is defined as 

x% = e ( r t p - o "7/ (*r - R t) > ( 3 °) 

ij 

where R^ xp and Rf 1 are the experimental values and theoretical predictions of the observables respectively and 
<j~ is the inverse of the covariance error matrix built from the statistical and systematic errors considering mutual 
correlations. 

We use 80 data points in our analysis; the total rate of the chlorine experiment (Homestake), the average rate of the 
gallium experiments (SAGE, GALLEX, GNO), 44 data points from the SK zenith-angle-spectrum and 34 data points 
from the SNO day-night-spectrum. The only correlation between the rates of the water Cerenkov experiments and 
radiochemical experiments is the uncertainty in the S B flux. Since we fit the shape of the spectrum for water Cerenkov 
experiments in our analysis, the covariance error matrix can be block diagonalized. We write contributions to x© 
from the rates of the radiochemical experiments, the SK zenith-angle-spectrum, and the SNO day-night-spectrum 
explicitly: 

X© — Xci.GaRatcs + X SK + XSNO (31) 

SK and SNO uses the pattern and intensity of the Cerenkov light generated by the recoiling electron in order to detect 
events due to electron scattering (ES): 

u x + e~ -» v x + e- (ES). (32) 

ES is sensitive to all neutrino flavors with reduced sensitivity to non-electron-neutrino components. Since SNO 
contains heavy water it is sensitive to charge-current (CC) and neutral-current (NC) reactions in addition to ES: 

v e + d -> p + p + e- , (CC), (33) 
v x + d -► n+p + v x , (NC). (34) 

The Cerenkov light generated by the recoiling electron is used to observe the CC events while the gamma ray from 
the neutron capture on deuterium is used to detect NC events. Time, location, direction, and energy of the CC events 
allow reconstruction of the solar neutrino spectrum. 

Since Cerenkov experiments have a higher threshold energy they are sensitive to only S B and hep neutrinos. The 
hep neutrino flux is much smaller than the S B neutrino flux, but hep neutrinos are somewhat more energetic. The 
production regions for these two components of the neutrino flux are not much different. For the sake of simplicity 
instead of dealing with two different sources of neutrinos, we add hep neutrino spectrum to the S B neutrino spectrum 

\ \ 4>B^b(E w ) + (j) hep \ hep (E v ) 

Ab{.& v ) ->■ t — — , (35) 



4>B -> 4>B + 4>hep, (36) 

and use them as a single source in the analysis of Cerenkov experiments. 

Cerenkov experiments are also live-time experiments. They can measure separate day and night rates or even divide 
their night rate into several zenith angle "a" bins. We may expect to see a different rate for each of these bins since 
MSW mechanism predicts earth regeneration effects at night when neutrinos pass through several layers of Earth 
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material. We incorporate those effects into our analysis by calculating the survival probability numerically at each 
zenith angle using a step function density approximation to the Preliminary Earth Model [H^. Survival probability 
for each zenith-angle bin is calculated by averaging the probability weighted with the exposure function "/(a)" of the 
detector. For SK we used exposure function given in Ref. 53]. For SNO we only used day-night bins and live-time 
information from Ref. |54|. For any zenith bin between a 3 min and oc J max , the "weighted" average survival probability 
(pS(y e ->u e ,E v ))ia: 



,E„)) = 



J°W f{a)P (v e -^v e ,E v ,a) da 
Ch* f(a)da 



(37) 



SK measures the kinetic energy of the recoiling electron and reports the data divided into several kinetic energy 
intervals. The kinetic energy assigned to the event by the detector is not always same as the true kinetic energy, 
instead it has a distribution around actual kinetic energy. This is characterized by a Gaussian shaped response 
function l55l: 



R(T, T 1 ) 

SSK 



1 



l\2 



—j= exp 

V27TS S K 

0.47^ MeV, 



(T - T") 



(38) 
(39) 



where ssk is the width of the Gaussian 56g , T 1 is the actual kinetic energy of the recoiling electron and T is the 
kinetic energy assigned to the same event by the detector (in Eqs. i|38|) and 139|) both T and T 1 are in MeV). By 
convolving differential ES cross sections of Ref. [5J7J with energy response function of the detector at each kinetic 
energy interval i, we get the "corrected" cross sections corresponding to the kinetic energy bin i: 



,{E V ) 



dT 



dT' R(T, T") 



da^(E„ T) 
dT 1 



(40) 



where T' max — E v /(1 + m e /2E v ) is the maximum kinetic energy that any electron can have due to kincmatical limits. 
After this step we have all the ingredients to calculate the rates for the zenith-spectrum bins, weighted survival 
probabilities for zenith-angle bins {{P 3 {v e — > v ei E v ))), and the corrected cross sections for kinetic energy spectrum 
bins {a\, v (E v )). If we consider only oscillations into active flavors the theoretical rate at SK for each zenith-spectrum 
bin is 



T>SK _ R SK 



L; i dE v X B (E v ) [ CT * c [E v ) {P' J [y e -> v e , E v , t)) 
- aUE v ){l-{pi{v e ^v e ,E v ,t)))], 



(41) 



where we use "fc" as a collective index for zenith-spectrum bins instead of the pair f. SK reports the ratio of 
number of observed events to the number of expected events under no oscillation condition. Dividing the above rate 
with SSM expected value (i.e. with the survival probability is 1) we obtain ratios to be compared with those given in 
Ref. 0. 

To calculate the error matrix for SK zenith-spectrum bins, experimental rates and uncertainties are taken from 
[58|. For each zenith and energy bin SK reports rates and statistical and systematic uncertainties. We take the 8 E> 
shape, energy scale and energy resolution uncertainties from which are actually calculated under no-oscillation 
condition and may result in discrepancies at higher confidence levels. An additional overall systematic offset error of 
2.75% is added to all bins. In calculating Xsk we introduce a free normalization parameter "ry" and minimize Xsk 
with respect to r\. In this way the total 8 B flux allowed to float freely. 



Xsk 



5> 

ij 

SNO response function for electrons is 

R(T, T') = 



uth\ —2 / r> exp 

T]R t )a lj (i? ■ 



rjR 



(42) 



1 



■ exp 



(T - T) 
2s 2 



27TSSNO 

ssno = (-0.0684 + 0.331 VT 7 + 0.0425T') MeV, 



(43) 
(44) 
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"Corrected" ES and CC cross sections can be obtained in a similar fashion as we did for SK (cf. Eq. f I40[l). NC 
events are mono-energetic. Neutrons are first thermalized and then captured on deuterium. All NC events originally 
have same energy, Tmc = 5.08 MeV. Their response function is 



R(T) = —= exp 

V27TSNC 



(T - T NC f 
2s 2 



(45) 



s NC = l.llMeV, (46) 
and the corresponding "corrected" cross sections are: 

a% c (E v ) = [™dTR{T) dG - c f- T) . (47) 

The response function of Eq. (|45|) spreads the neutral current events in energy. We use the "forward-fitting" technique 
described in [54| to calculate x"sno- Event rates for each type of reaction are 

^f 5 = Rfj = <t>B J dE v \ B {E v )[ai B {E v ){Pi e {E v ,t)) + ai^E v ){l - (P^E^t)))], (48) 

Rk° = Rff = 4>b j dE u X B (E U ) [ah c (E u )(Pi e (E u , t))] , (49) 

R% c = Rff = <j> B f dE v X B {E U ) [<r* N c(Ey)] • (50) 

In our calculations we used the neutrino-deuteron cross-sections of Ref. [5^ calculated using the effective field theory 
approach. We fixed the counter-term, L\a, of this approach so that it reproduces the calculation of Ref. which 
incorporates the first-forbidden matrix elements in the calculation of the neutrino-deuteron cross sections. Theoreti- 
cally expected rate is calculated by adding background contributions like the so-called Low Energy Background (LB) 
and Neutron Background (NB) to the sum of CC, NC and ES events. We take these background contributions from 
Ref. |5J|. The expected event rate for each bin is 

nth _ r/ES i jyCC , nJVC , r>LB i r>NB /n \ 

In calculating error matrix for SNO zenith-spectrum bins, experimental rates and uncertainties are taken from Ref. 
[54|. Statistical errors are calculated from the data reported by SNO and systematic uncertainties (shape scale and 
resolution) are taken from Ref. |49j |. Other systematics like vertex accuracy, neutron-capture efficiency, etc. are taken 
from Ref. In calculating Xsno we multiply sum of CC, NC and ES events (without backgrounds) by a free 

normalization parameter "77" and minimize Xsno with respect to r\ as we did for SK (cf. Eq. I|42|l . In this manner 
total 8 B flux allowed to float freely without affecting backgrounds. 

The last component of x% is from the radiochemical experiments. In the evaluation of error matrix for Xc\ Ga Rates 
we follow the procedure described in Ref. [6l| . 

In the presence of oscillations, energy averaged cross section of neutrinos from source m at detector i can be 
calculated by convolving the neutrino spectrum X m (E u ) for the corresponding source, the cross section <Ji{E u ) for the 
corresponding detector and the survival probability {P{v e — > v e ,E v )) (averaged over source distributions in the Sun 
and weighted with the exposure function of each detector) 

C^i = J dE v \ m {E v ) a t {E v ) (P(u e -» u e , E v )). (52) 

Then event rate at detector i is simply: 

R l h = R ™i = ^ mC ™ ( 53 ) 

m rn 

SSM neutrino fluxes <p m depend on the SSM input parameters X/~. The correlations between neutrino fluxes are 
parameterized by the logarithmic derivative: 

dln<bf SM 

a * - ~dhTx^~ (54) 
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AlnXfe and AlnCfci are la relative errors of SSM input parameters and energy averaged cross sections respectively. 
We adopt the values of these parameters from Ref. |61j . 
From all above, one can write 

4 = 5 tj o^o^ + Si, Ru r K r ( A ln 2 + Y, R % Rt £ J2 amkank ( A lnX ^) 2 • 

k ran k 

For global analysis we take: 

2 _ 2 , 2 

Xdobal — X © T XKamLAND 

KamLAND detects reactor neutrinos in 1 kiloton of liquid scintillator through the reaction: 



(55) 



(56) 



P + Ve 



n + e 



+ 



(57) 



We use the phenomenological parameterization of the energy spectrum of the incoming antineutrinos given in Ref. 
Ill: 



ao+aiE v +a 2 E, y 



(58) 



where the parameters vary with different isotopes. The spectra of antineutrinos coming from each detector, cf>i(t, E u ), 
can be calculated using the thermal power and the iso trop ic composition of each detector. The effects of incomplete 
knowledge of the fuel composition are explored in Ref. 63] . We use the time-averaged fuel composition for the nuclear 
reactors given by the KamLAND collaboration [llj . 

Survival probability for electron antineutrinos coming from the reactor i is: 



sin #i3 + cos 6*i3 



i • 2 ozi -2 ( l-27Am 21 i,- 

1 — > sin 2#i2 sin 1 



E v 



(59) 



where Li are the reactor-detector distances. We denote the energy resolution function of KamLAND by R(TT') 
where T, T 1 are the observed and the true positron energies. The energy resolution is given as 7.5%/y/E(MeV) 
The number of expected events for each bin at KamLAND can be calculated by convolving the cross section, <j(E u ) : 
with the reactor spectra, survival probabilities and the resolution function of KamLAND: 



A* h = J dE v a{E v ) ^j^- Pjfre - E v ) J dTR{T, T') . 



(60) 



where the electron kinetic energy is T' ~ E v — 0.8 MeV, and a{E v ) is the lowest order cross section given in Refs. 

His!: 



<r{E v ) 



2tt 2 



m\jr n 



p e E e 



(61) 



in which E e = T' + m e , f = 1.69 is the integrated Fermi function for neutron beta decay, and r„ is the neutron 
lifetime. 

KamLAND reports its results in 13 bins above the threshold. Due to low statistics, we use the prescription of [6|J 
in analyzing KamLAND data. We write 



X 2 KamL AND -Sped ~ 2 ^ 



(ryA^ - AT. p ) + N? 



'In 



r,Nl 



Hi 



(62) 



sys 



where <J 2 sys = 6.75% is the systematic uncertainty, and minimize the sum in Eq. 116211 with respect to r\. For the total 
rate analysis we use 



X KamLAND— Rate 



(_/V t/l — _/v ea:p ) 2 



(63) 



with errors added in the quadrature to calculate <7j|. N exp is given in 11] and N th is calculated at each value of 
oscillation parameters similar to the binned expected event rates. 
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IV. RESULTS AND CONCLUSIONS 

In our calculations we use the neutrino spectra given by the Standard Solar Model of Bahcall and collaborators [j| . 
It it is numerically more convenient to follow the evolution of the matter eigenstates in the Sun and in the Earth (or 
mass eigenstates in vacuum), a procedure which we adopted. We take into account the distribution of various neutrino 
sources in the core of the Sun and resulting non-linear paths of the neutrinos. Thus neutrinos comin g fr om the other 
side of the Sun may have double resonances. In the Sun we used the Landau-Zener approximation [H. 1671 168], We 
divided the Sun into several shells which were in turn divided into several angular bins, calculated the derivative of 
the electron density in the Landau-Zener approximation numerically for both radial and non-radial neutrino paths 
and averaged the survival probabilities over the initial source distributions. 

Survival probabilities in the Earth depend on zenith angles. As was described in the previous Section we adopted 
the Preliminary Earth Model density 52] to solve neutrino evolution oscillations numerically in the Earth. 

In our calculations we ignore the possibility of density fluctuations in the Sun l69l F(i| . Recent data indicate that 
such fluctuations are less than one percent of the Standard Solar Model density [7l|. We similarl y ignore p ossible 
mixing of sterile components. Sterile neutrinos can play a very important role in supernova r-process |72ll73ll74| or the 
big-bang nucleosynthesis [T^. It is worth emphasizing that active-sterile mixin g ca n be too small to be detectable in 
solar neutrino experiments (for a discussion of various possibilities see e.g. Ref. |7(j) and yet may still have significant 
astrophysical impact. 

We first present calculations where we took the value of 6*13 to be zero and considered only the solar neutrino 
data. Allowed regions of neutrino parameter space when each solar neutrino experiment is considered separately are 
shown in Fig. [2 One observes that either Sudbury Neutrino Observatory or SuperKamiokande individually already 
significantly limit the neutrino parameter space. Our results agree well with the own analyses of these experimental 
groups. Allowed regions of neutrino parameter space when all solar neutrino experiments are combined together are 
shown in Fig. [21 The combined solar neutrino data already rule out the so-called LOW region where 5m\ 2 < 10~ 7 eV 2 
at the 99 % confidence level. We find the best fit (minimum x 2 ) values of neutrino parameters to be tan 2 812 ~ 0.45 
and 5ra\ 2 ~ 7.08 x 10 -5 eV 2 from our combined analysis of all the solar neutrino data. Our minimum \ 2 value is 
67.2 for 80 data points and 2 parameters. 

We next turn our attention to KamLAND data while still keeping the value of #13 to be zero. In the upper panels of 
Fig. 0]we display allowed regions of the neutrino parameter space from the KamLAND data only. (Results with using 
the total rate only is at the left-hand side and results with using the binned data, which provide more information 
about the neutrino spectrum, is at the right-hand side). The result of our global analysis combining the solar date 
with data from KamLAND is shown at the lower right-hand side panel. For convenience of presentation we do not 
show the lower values of Sm 2 2 on the graph, but the LOW solution is completely eliminated. For comparison we 
re-plot the neutrino parameter space obtained from the solar neutrino data only at the lower left-hand side of the 
Figure. KamLAND data significantly shrinks the LMA region (the lower right-hand side of the Figure). We find 
the best fit (minimum x 2 ) values of neutrino parameters to be tan 2 6*12 ~ 0.46 and Sm 2 2 ~ 7.1 x 10~ 5 eV 2 from 
our combined analysis of all the solar neutrino and KamLAND data. Our minimum \ 2 value is 73.2 for 80 + 13 
data points and 2 parameters. Inclusion of the KamLAND data does not noticeably change the best fit values of 
the neutrino parameters, however KamLAND, being a terrestrial experiment with very different statistical errors, 
provides a completely independent test of the results from the solar neutrino experiments. Fig. ^also illustrates that 
mixing of the (solar) neutrinos and (reactor) antineutrinos are very similar, very likely to be identical. 

It is instructive to investigate how well the neutrino parameter space is constrained. To this extend we plot allowed 
regions of the neutrino parameter spaces obtained by combining data from a single solar neutrino experiment with the 
KamLAND data in Fig. In calculating the parameter space shown in this Figure we continued to take the value 
of 6*i3 to be zero. We observe that any single solar neutrino experiment taken together with KamLAND significantly 
constraints the parameter space. The LOW region, which is not shown on these plots, is again completely eliminated 
for each case. Real-time Cerenkov detectors are slightly more constraining than the radiochemical experiments in this 
regard. It is also interesting to realize that one no longer needs all the solar neutrino experiments to determine the 
neutrino parameters. We are getting closer to realizing the initial goal of the solar neutrino experiments, eloquently 
stated in the seminal papers of Bahcall and Davis |77j , namely to use solar neutrino data to better understand the 
Sun. (For a preliminary effort see Ref. |45jl- 

We next examine the effects of a non-zero value of $13. In Fig. 0we show how the parameter space changes as a 
function of #13 when 812 is kept fixed. Here we use Eq. (|29|l to calculate the 3-flavor neutrino survival probability 
and perform the \ 2 analysis for three parameters (#12, #13, and 5m\ 2 ). We find the best fit (minimum % 2 ) values 
of neutrino parameters to be tan 2 812 ~ 0.46, cos 4 6*13 ~ 1, and 8m\ 2 ~ 7.1 x 10 -5 eV 2 . Note that the confidence 
level regions in the pair that corresponds to the best fit (the lowest right-hand side pair in Fig. [BJ are larger than 
the corresponding pair obtained with 2-flavor analysis (lower pair in Fig. @J since as one goes from the former to the 
latter the number of parameters are reduced by one. To illustrate the change in the quality of the fit as tan 2 6*12 is 
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changed we present confidence levels for a several values of tan 2 #12, some of which are clearly very far away from the 
optimal solution. It is interesting to note that the combined Solar neutrino and KamLAND date provide a limit of 
cos 4 8i3 > 0.8 at the 90 % confidence level. We are able to put such a limit since we investigated confidence levels 
for larger values of #13, otherwise solar data alone are not sufficient to constraint this angle (cf. Ref l78T0. T his limit 
currently is not as good as the one obtained from the completed reactor disappearance experiments |30L l31| ; however 
the situation may change after a few years of data taking at KamLAND [xjj . 
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FIG. 1: Several components of the electron neutrino survival probability given in Eq. J1IJ. In this figure 6 13 = tt/6, 0i 3 = tt/20, 
E v = 1 MeV, and (Jm^ = 10 -3 eV 2 . The upper solid line is the total survival probability. The filled squares on top of the 
figure represent the contribution of the term proportional to cos 4 #13 alone to the survival probability. The middle long-dashed 
line is the contribution of the term proportional to sin 4 #13 alone to the survival probability. The lower dash-dotted line is the 
contribution of the correction terms proportional to £ only. 
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FIG. 2: Allowed regions of the neutrino parameter space when each solar neutrino experiment is considered separately. In this 
figure #13 is taken to be zero. The shaded areas are the 90 % confidence level regions. 95 % (solid line), 99 % (log-dashed line), 
and 99.73 % (dotted-line) confidence levels are also shown. 
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FIG. 3: Allowed regions of the neutrino parameter space when all solar neutrino experiments (chlorine, all three gallium, SNO 
and SK experiments) are included in the analysis. In this figure #13 is taken to be zero. The shaded area is the 90 % confidence 
level region. 95 % (solid line), 99 % (log-dashed line), and 99.73 % (dotted- line) confidence levels are also shown. 




FIG. 4: Allowed regions of the neutrino parameter space with #13 taken to be zero. In the upper-left panel only KamLAND 
total rate and in the upper-right-hand panel binned KamLAND data are used. Lower-left hand panel depicts the LMA solution 
obtained from only the solar neutrino experiments. All solar neutrino experiments are combined with the KamLAND date to 
obtain the lower right-hand panel. The shaded areas are the 90 % confidence level regions. 95 % (solid line), 99 % (log-dashed 
line), and 99.73 % (dotted-line) confidence levels are also shown. 
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FIG. 5: Allowed regions of the neutrino parameter space obtained when data from a single solar neutrino experiment is combined 
with the KamLAND data. The shaded areas are the 90 % confidence level regions. 95 % (solid line), 99 % (log-dashed line), 
and 99.73 % (dotted-line) confidence levels are also shown. 
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FIG. 6: Allowed regions of neutrino parameter space when both mixing angles O12 and 813 are varied. Six pairs of Sm^ 
vs. cos 4 S13 plots are presented for representative values of tan 2 612 (0.3, 0.4, 0.43, 0.46, 0.5, 0.7, and 0.8). The figures on 
the left-hand side of each pair are obtained from the combination of all solar neutrino experiments while the figures on the 
right-hand side of each pair also include the KamLAND data. The shaded areas are the 90 % confidence level regions. 95 % 
(solid line), 99 % (log-dashed line), and 99.73 % (dotted-line) confidence levels are also shown. 



